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Abstract 

We construct an analytical model derived from nuclear reaction theory and hav- 
ing a simple functional form to demonstrate the quantitative agreement with the 
measured cross sections for neutron induced reactions. The neutron-nucleus total, 
reaction and scattering cross sections, for energies ranging from 5 to 700 MeV and 
for several nuclei spanning a wide mass range are estimated. Systematics of neu- 
tron scattering cross sections on various materials for neutron energies upto several 
hundred MeV are important for ADSS applications. The reaction cross sections of 
neutrons are useful for determining the neutron induced fission yields in actinides 
and pre-actinides. The present model based on nuclear reaction theory provides 
good estimates of the total cross section for neutron induced reaction. 

PACS numbers: 24.10.-i, 25.40.-h, 28.20.-v, 28.20. Cz, 28.20. Fc, 24.10. Ht, 25.60.Dz 
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1 Introduction 



The cross sections for neutron induced reactions by nuclei at very high en- 
ergies upto several hundred MeV are required in a number of fields of study 
in basic science as well as many of applied nature such as radioactive ion 
beam (RIB) (Diamond 1999; Essabaa 2003) production or the accelerator- 
driven sub-critical systems (ADSS) and their neutronics (Boudard et al. 2000; 
Broeders and Broeders 2000; Demirkol et al. 2004; Kaplan et al. 2009; Yapycy 
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et al. 2007; Demirkol et al. 2008). The sub-critical reactor is driven critical by 
spallation neutrons produced by bombarding high energy proton beam with 
high current (greater than ~ 10mA) (Rubbia et al. 1995) on a heavy element 
target. Such a system serves a dual purpose of energy multiplication and waste 
incineration of the long lived radioactive waste produced in reactors based on 
thermal neutron induced fission. Unlike the thermal neutron induced fission, 
the energy spectrum of spallation neutrons can reach up to several hundred 
MeV. In this context it is important to study the systematics of neutron ab- 
sorption and scattering cross sections on various nuclei for neutron energies 
up to several hundred MeV. 

Often, these cross sections are evaluated using phenomenological optical po- 
tentials and much effort has gone into defining global sets of parameter values 
for those optical potentials with which to estimate cross sections as yet unmea- 
sured. In a recent study, Koning and Delaroche (Koning and Delaroche 2003) 
gave a detailed specification. In this context it would be utilitarian to study 
the systematics of neutron absorption and scattering cross sections on vari- 
ous nuclei well approximated by a simple convenient functional form (Deb and 
Amos 2004; Deb, Amos, and Karataghdis 2004). Indeed that is so and we show 
herein that the simple analytical model can be used to estimate these cross 
sections without recourse to optical potential calculations which are limited 
up to 150-200 MeV. 



2 Theoretical formalism 

Although the optical model calculations give reasonable account of all the 
data presented here, these calculations involve many partial waves and it is 
not easy to get a simple interpretation or intuitive picture of the basic processes 
involved. The basic picture of the nuclear Ramsauer model is that the forward 
scattering amplitude for a neutron incident on a nucleus is given by 



where the complex phase shift 5i may be considered as independent of angular 
momenta I of the Ith partial wave (Lawson 1953), so that the above expression, 
after summing over / upto a maximum value L = kRch (Angeli and Csikai 
1970; Angeli and Csikai 1971; Angeli, Csikai and Nagy 1974), reduces to 
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where A — 1/k — h/\/2mE, E is the incident neutron energy in the center 
of mass system, m is the reduced mass of the neutron-nucleus system and 
Rch is the channel radius beyond which partial waves do not contribute. The 
quantity e^*^' is replaced by ae*'^ with (3 being two times the real part of the 
phase shift and a accounting for the attenuation or loss of flux (arising out 
of imaginary part of the phase shift) . From partial wave analysis of scattering 
theory, we know the standard expressions for scattering and reaction (T^ 
cross sections as 

asc = ^^i {21 + 1)|1 - Vll'': <^r = ^^i {21 + 1)[1 - Ivif] (3) 



where the quantity rji = e With the assumption that the phase shift 6i is 
independent of / and the summation over partial waves / is upto kRch only, it 
follows that 



asc = 7r(i?ch + A)^(l + q;^ - 2q;cos^), 
ar^n{Rch + A)^{l-a^), 

CTtot = CTsc + Crr^ 2TT{Rch + A)^(l " « COS /3) (4) 

where (3 = 2Re5; = 2ReS, a = e"^^™''' = e"^^™'' and summing over I from to 
kR,h yields {21 + 1) = {kR,h + if = k\R,h + A)'- 

The optical potential for nuclear n-N interaction can be written as —V — iW 
with V and W as positive quantities, and contains no Coulomb interaction. 
The phase shift 5 in a WKB approximation is [/ K'dr — J k'dr] and the real 
part of it to a zeroth order approximation (Mohr 1957) for a square well with 
radius R is {K — k)R where K is the real part of K' and k' = k is real due 
to absence of potential. The real wave numbers inside and outside the nucleus 
are, therefore, given by K'^ — 2m{E + V)/Ti^ and k"^ — 2mE/U^ respectively. 
Hence ^ is determined by the real potential V, 

I3 = 2{K -k)R = 2^^^[^/ETV -^/E]R (5) 



whereas the attenuation factor a is determined primarily by the imaginary 
potential W, 
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where A is the mean free path of the neutron inside the nucleus. The average 
chord length ^ of a neutron passing through a nucleus can be derived as 

R ^ A 

j 2Vi?2 _ x'^{I2T:xdx) I J I27rxdx = -R (7) 



where / is the neutron flux that is the number of neutrons incident per unit 
area. Since R (x As (i? ~ ro^s), the above arguments imply that 

P = [Ve + V - Ve] (8) 



where Pq = ^^"'•^"^•'^ whose value is approximately 0.6, and the attenuation fac- 
tor which is much less than unity but increases with energy (which is obvious 
from its expression) is given by 



a 



exp[-aoroA^W/VE + V] (9) 



where ckq = ^ whose value turns out to be 0.2929 and Tq is the nuclear 
radius parameter. The first term Va = Vq + Vi(l — 2Z/A) + V2/A of the 
real potential V = Va + VeVE contains both the isoscalar and the isovector 
(Lane 1962; Satchler 1983) components of the optical potential (Gould 1986; 
Anderson and Grimes 1990) where Z is the atomic number of the target 
nucleus, whereas the second term accounts for its energy dependence. The 
imaginary potential W is taken as = Wq + WeVE + V since the total 
kinetic energy of the neutron inside the nucleus with attractive potential well 
of depth V is E + V. As the magnitude of the real part of the optical potential 
decreases with energy while the same for imaginary part increases, this implies 
that Ve is negative whereas the We is positive. 

It is worthwhile to mention that in all the previously published papers, the 
expression for the phase shift (3 has been explained as caused by refraction 
in a sphere of radius R with refractive index n like Peterson (Peterson 1962). 
We point out here that the derivation of ref. (Peterson 1962) is not profound 
and provided above an alternative, theoretically concrete, derivation for the 
phase shift using WKB method of quantum mechanics as applied to scattering 
by a spherical potential well. The drawback of the derivation of ref. (Peterson 
1962) is that the neutron (although massive) is treated like a photon and as 
its velocity inside nucleus and vacuum are proportional to y/E + V and V^, 
respectively, it would result in bending of the ray (as in optics) away from 
the normal inside nucleus (in fig. 16 of ref . (Peterson 1962) it is shown just 
the opposite) where velocity is more. This would lead to the existence of the 
critical angle sin~^JE/{E + V) beyond which there is no transmission (even 
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in an attractive nuclear potential) and a refractive index less than vacuum 
for the nuclear medium which are physically unacceptable. Even then, if one 
sticks to Peterson's assumption of a light ray, the average chord length inside 
nucleus, with ray bending away from the normal, turns out to be less than 
our result of 4i?/3 as opposed to greater than 4i?/3 as derived in ref. (Peterson 
1962). However, these results go over to our result of 4/?/3 for refractive index 
equal to one, that is, reach our result asymptotically at energies higher than 
magnitude of the real part of the nuclear potential. 



3 The analytical model calculations 

The Ramsauer model can be fitted to the experimental neutron total cross 
sections using Eq.(4). The radius of the nuclear potential is given hy R = roAs 
whereas the channel radius can be parametrized (Gowda and Ganesan 2006) 
as Rch — roAs +rAy/E + r2 with = rio InA + rn/hiA. The Ramsauer model 
fits yield rio = -22.98 x 10-^ rn = 10.27x IQ-^ r2 = 23.22x10-2, I/q = 46.51, 
Vi = 6.74, V2 = -117.52, Ve = -3.22 and /3o = 0.5928. These values are very 
close to or within the limit of the parameter values obtained in our earlier work 
(Mukhopadhyay, Lahiri and Basu 2010). The value of cto is kept fixed at 0.2929 
and the non linear least square fits yield the values for the imaginary potential 
Wo = 5.293 MeV and its energy dependence We = 33.88 x 10"^. The nuclear 
radius parameter Tq is also fitted reasonably well to 1.378 A'^ fm which means 
that the nuclear potential radius R = r'^A^'^'^ where 7 = 7.93 x 10^'^ is a very 
small number (needed for fine tuning) compared to ^. These extracted model 
parameters provide global fits to the neutron total cross sections spanning 
quite a large number of target nuclei. In our earlier work (Mukhopadhyay, 
Lahiri and Basu 2010) the effect of the imaginary potential was included in 
the parameter a in somewhat ad-hoc manner lacking justification of it to 
be weakly mass dependent. The present work is an improvement over our 
earlier work where the imaginary part of the nuclear potential is now treated 
appropriately with explicit appearance of the imaginary part of the nuclear 
potential. 



4 Results 

Each calculation is performed at neutron incident energy intervals of 1 MeV 

and for various elements. In Figs. 1-15, the variations of the total cross section 
{cTtot), the scattering cross section (c^c) and the reaction cross section (o-^.) with 
incident neutron energy are plotted for ^^^U, ^^^Th, ^"^Bi, ^"^Pb, ^^^Au, ^^^W, 
184^, 182^^ isirp^^ 93^^^ 902r, 59^^^ ss^n, ^o^a and ^^^P target nuclei. The 
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continuous lines represent the total cross sections, the dashed lines represent 
the scattering cross sections and the dotted lines represent the reaction cross 
sections while the hollow circles represent the experimental data (Finlay et 
al. 1993; Abfalterer et al. 2001; Dietrich et al. 2003) for atot which show good 
agreement with the measured cross sections. In Fig. 16 estimates of these cross 
sections are plotted for ^^^Pu target. 



5 Summciry and conclusion 

In summary, we constructed an analytical model, justified it from the opti- 
cal model and nuclear reaction theory approach and applied it to derive the 
systcmatics and performed calculations of neutron-nucleus total cross section 
(o-tot), the scattering cross section (asc) and the reaction cross section ((7^) and 
predicted these cross sections for heavier actinides. The extracted parameters 
for the present analytical model provide global fits to the neutron total cross 
sections spanning quite a large number of nuclei. We conclude that the present 
estimates of neutron scattering cross sections are very important for the re- 
actor physics calculations Diamond 1999; Essabaa 2003) production or the 
accelerator-driven sub-critical systems (ADSS) and their neutronics (Boudard 
et al. 2000; Broeders and Broeders 2000; Demirkol et al. 2004; Kaplan et al. 
2009; Yapycy et al. 2007; Demirkol et al. 2008; Rubbia et al. 1995) of the 
ADSS applications. Moreover, the neutron-nucleus reaction cross sections are 
very useful for the theoretical calculations of Radioactive Ion Beam (Dia- 
mond 1999; Essabaa 2003) production. These can also be used for performing 
Hauser-Feshbach (Hauser and Feshbach 1952) calculations with Monte-Carlo 
simulations (Pace2 code 1984) to estimate the cross sections for neutron in- 
duced fission, evaporation residues or evaporation neutron multiplicities and 
for comparison of the photon (Mukhopadhyay and Basu 2007; Mukhopadhyay 
and Basu 2009) versus neutron induced fission as well and may serve as inputs 
to intranuclear cascade codes such as MCNPX package. 
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Figure captions 

Fig.l. The plots of total cross section {atot), the scattering cross section (a^c) 
and the reaction cross section ((7^) versus incident neutron energy for ^^^U 
target. The continuous line represents the total cross section, the dashed line 
represents the scattering cross section, the dotted line represents the reaction 
cross section and the hollow circles represent the experimental data. 
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Fig. 16. The plots of total cross section {(Jtot), the scattering cross section 
(cTsc) and the reaction cross section (cXr) versus incident neutron energy for 
^^^Pu target. The continuous hne represents the total cross section, the dashed 
line represents the scattering cross section and the dotted line represents the 
reaction cross section. 
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